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Abstract: We consider the minimal seesaw model in which two gauge singlet right handed 
neutrinos with opposite lepton numbers are added to the Standard Model. In this model, 
the smallness of the neutrino mass is explained by the tiny lepton number violating coupling 
between one of the singlets with the standard left-handed neutrinos. This allows one to have 
the right handed neutrino mass at the TeV scale as well as appreciable mixing between 
the light and heavy states. This model is fully reconstructible in terms of the neutrino 
oscillation parameters apart from the overall coupling strengths. We show that the overall 
coupling strength for the Dirac type coupling between the left handed neutrino and one 
of the singlets can be restricted by consideration of the stability bounds on the electroweak 
vacuum. Incorporating this bound, the overall coupling strength of the small lepton number 
violating coupling can also be constrained from neutrino oscillation data. In this scenario 
the lepton flavor violating decays of charged leptons can be appreciable which can put 
further constraint on yy, for right-handed neutrinos at TeV scale. We discuss the combined 
constraints on y^ for this scenario from the process ^ — >• 67 and from the consideration 
of vacuum stability constraints on the Higgs self coupling. We also briefly discuss the 
implications for neutrinoless double beta decay and possible signatures of the model that 
can be expected at colliders. 
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1. Introduction 

This year marks the announcement of a very significant result in the field of particle physics. 
The ATLAS and CMS collaboration of LHC experiment has reported the observation of a 
new neutral boson. The mass of this particle is reported to be [1,2] 

M = 126.0 ± 0.4 ± 0.4 GeV (ATLAS); M = 125.3 ± 0.4 ± 0.5 GeV (CMS) (1.1) 

This, in all probability, is the Standard Model Higgs boson which has eluded scientists so 
far. Further analysis and data would confirm this and would also explore if there is any 
hint of new physics beyond the Standard Model. 

The Higgs boson is responsible for giving mass to all the fundamental particles. How- 
ever how neutrinos get their mass still remains an enigma. The existence of neutrino masses 
and flavor mixing are already established by the oscillation experiments. The oscillation 
of three known neutrinos are characterized by two mass squared differences with values 
~ 10~^ eV^ and ~ 10~^ eV^. There is also a mass bound on the sum-total of the light 
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neutrino masses from cosmology [3-5], which is < 0.5 cV. The mass squared differ- 

ences inferred from oscillation data along with the cosmological mass bound indicate that 
the neutrino masses are much smaller than the corresponding charged fermion masses. 

The most natural mechanism which can explain such small masses is the Seesaw Mech- 
anism which postulates a heavy particle at some high scale determined by the mass of this 
particle. The tree level exchange of such a particle gives rise to an effective neutrino mass 
operator kIlIl(/)(/)/M in terms of the Standard Model lepton and Higgs doublet fields [6]. 
For a typical Higgs vacuum expectation value of ~ 100 GeV and a typical value of the 
heavy scale M ~ 10^^ GeV the above operator implies a neutrino mass of 0.1 eV. Since 
this term violates lepton number by two units the smallness of the neutrino masses gets 
correlated with the Majorana nature of the neutrinos. The most popular way to have an 
ultraviolet complete theory to generate the above dimension 5 effective operator consists 
of adding a heavy singlet state to the Standard Model which corresponds to the so called 
type-I seesaw mechanism [7-10]. This is already a step beyond the Standard Model of par- 
ticle physics. Since the LHC started operation, a natural question which has been explored 
in the literature quite extensively is the possibility of observing signature of seesaw at the 
LHC. This will require the mass of the heavy particles to be of the order of TcV scale. 
For the canonical type-I seesaw mechanism arising from SO(IO) theories, this is difficult 
and one has to appeal to cancellations coming from flavor symmetries [11-13]. Another 
option to relate small neutrino mass to TeV scale physics is the inverse seesaw mechanism 
which was first considered in [14] . In models with inverse seesaw mechanism, one includes 
additional singlet states with lepton number opposite to that of the right handed neutri- 
nos. The Majorana mass of the singlet breaks lepton number softly. The smallness of the 
neutrino mass can be related to the smallness of this parameter. For the canonical seesaw 
the neutrino mass is inversely proportional to the mass of the heavy particle which deter- 
mines the scale at which lepton number is violated. While in this case the neutrino mass is 
proportional to the mass of the singlet characterizing the lepton number violation. Hence 
this is termed as inverse seesaw. This is considered to be more natural because as this pa- 
rameter tends to zero the neutrino mass also becomes zero and lepton number symmetry is 
reinstated. This class of models with inverse seesaw mechanism is testable in various ways. 
Non-unitary mixing between light and heavy particles can be large and can be probed at 
colliders [15-17]. Future neutrino factories may also be sensitive to such mixings [18-20]. 
Lepton flavor violating processes can be appreciable [21,22]. The non-unitary effect can 
also play a non-trivial role in relating CP violation responsible for leptogenesis with low 
energy CP violation [23,24]. 

One variant of the inverse seesaw models of the type discussed in [14] is the linear seesaw 
model [25-27]. In this class of models, a tiny lepton number violation is incorporated 
through the term that couples the left-handed neutrinos and one of the singlet states. 
Within this scheme the minimal model consists of three left-handed neutrinos and just 
two singlet states. This is similar to the minimal type-I seesaw model [28-31], but in 
this case the two singlet states are assigned opposite lepton numbers. Thus only one 
of the singlets has a lepton number conserving Dirac type coupling with the left-handed 
states. The coupling of the other singlet with the left-handed states violate lepton number 
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symmetry and in the limit this couphng is zero the Lagrangian acquires an extra global 
U{1) symmetry. Therefore this should be naturally small. The smallness of the neutrino 
masses can now be related to the smallness of this coupling and the mass of the singlets 
can be at TeV leading to phenomenological consequences as in the inverse seesaw model. 
The advantage in this case is, the minimality of the scheme allows the reconstruction of 
the Yukawa matrices in terms of oscillation parameters apart from the overall coupling 
strengths [32]. 

However, assuming the particle observed in CMS and ATLAS is the Higgs Boson 
and its mass to be as given in Eq. (|l.l| ) opens up an avenue for constraining the Dirac 
Yukawa couplings in seesaw models from the consideration of stability of the electroweak 
vacuum [33,34]. It is well known that because of quantum corrections, the Higgs self- 
coupling A diverges for higher values of Higgs mass and goes to negative for low values of 
Higgs mass near Planck-scale (Mpi). Assuming no new physics between SM and the Planck 
scale, Higgs mass was found to be in the range 126 — 171 GeV for A(at Mpi) to be in the 
range [0,7r] [35,36]. The upper bound called the "triviality bound" essentially embodies 
the perturbativity of the theory. The lower bound known as the "vacuum stability bound" 
is obtained from the fact that a negative A makes the potential unbounded from below 
and the vacuum would be unstable. The presence of new Yukawa couplings in seesaw 
models modifies the /3 function of Higgs self-coupling. In the conventional type-I seesaw 
model, generation of small neutrino mass requires the mass scale of the singlet to be of 
the order of 10^^ GeV for Dirac Yukawa Coupling Yi, ~ C(l)- It was observed in [33] 
that the presence of this extra coupling increases the lower bound of the Higgs mass from 
vacuum stability constraints, gradually reaching the perturbativity bound. However as the 
mass scale of the heavy field is lowered, Yi, has to become less in order to get nii, ~ 0.1 
eV and below a certain value of the mass scale the additional contribution does not play 
any significant role. However from the point of view of relevance at LHC many models 
have been considered in the literature which can give rise to small neutrino masses with a 
relatively large Yukawa coupling even with the heavy field at the TeV scale. Hence in such 
models the effect of the Yukawa term can be significant in the running of A. Moreover, as 
the neutrino Yukawa runs from TeV to Planck scale, the effect can be large [37,38]. Since 
the presence of this term drives A towards a more negative value, the imposition of the 
condition \(Mpi) = 0, enables one to put constraint on the Yukawa coupling strength if 
the Higgs mass is known [33,34,37,38]. 

In this work we consider the Minimal Linear Seesaw Model (MLSM) which can natu- 
rally accommodate TeV scale singlets. We show that it is possible to constrain the unknown 
Dirac- Yukawa coupling strength y^/ in this model from the considerations of vacuum sta- 
bility of the scalar potential. Recently, vacuum stability bounds on the Dirac Yukawa 
coupling in TeV scale seesaw model have been obtained in [37,38]. However, in canoni- 
cal seesaw models one needs to make some assumptions about the structure of the Dirac 
type Yukawa matrix Y^, and the right handed Majorana mass matrix M^. On the other 
hand, for MLSM this is already completely determined in terms of oscillation parameters 
apart from the overall Yukawa coupling strengths and hence one need not make any further 
assumptions on the structure of the mass matrices. This feature makes it particularly suit- 
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able for studying vacuum stability constraints. From the bounds on i/i, thus obtained, the 
small coupling t/s can be constrained using the measured values of mass squared differences 
in neutrino oscillation experiments. Thus the overall coupling strengths which were the 
undetermined parameters in MLSM now gets constrained. Since the heavy singlet states 
in this model are at TeV scale, lepton flavor violating decays of charged leptons are not 
suppressed and from the bound on the branching ratios of these processes it is possible to 
constrain Uu/Mr. We consider the bound on the process ji ^ e-^ and discuss the upper 
bound obtained on together with the constraints from vacuum stability as a function of 
the mass scale Mr. We also comment on the implications of this model for neutrinoless 
double beta decay and discuss the possible collider signatures. 

The plan of the paper is as follows. In the next section we discuss the minimal singlet 
seesaw mass matrix. Section 3 describes the running of the self coupling A and investigates 
if the stability of the potential upto the Planck scale can be achieved assuming Higgs mass 
in the range measured by ATLAS and CMS experiment. How this condition changes in the 
minimal linear seesaw model is also examined. We also obtain the constraints on from 
vacuum stability in MLSM for known values of Higgs masses. In the subsequent section 
we delineate the bounds on the coupling yg from the vacuum stability bound on y^ and 
neutrino oscillation data. In section 6 we study the phenomenological implications of this 
model for charged Lepton Flavor Violation (LFV). In section 7 we consider neutrinoless 
double beta decay (Oi^/3/3) in this model. In the next section we comment on possible 
collider signatures and finally present the conclusions in section 9. 



2. Singlet Seesaw Models 

The Yukawa part of the most general Lagrangian involving extra singlet states can be 
written as 

-£ = NrY,4)Hl + SYs4>Ul + 'SMrN% + ^^S"" + ^iV^M^iV^ + h.c. 

where II = (z^xj^c)^, x = e,iJ,,T. Il,Nr and S have lepton number 1,1,-1, respectively. 
After spontaneous symmetry breaking, the (p field acquires a vacuum expectation value 
{v/^/2) and Yi,v/y/2 = mn gives rise to the Dirac mass term while the term Ysv/V2 = ms 
breaks lepton number. In the above Lagrangian lepton number violation stems from the 
terms with coefficients Yg, jj, and Mjv and thus the symmetry of the Lagrangian is enhanced 
(lepton number becomes an exact symmetry) in the absence of these terms. Therefore these 
co-efficients are expected to be naturally small (i.e. there is no fine tuning or unnaturalness 
in keeping these terms to be very small) according to 't Hooft's naturalness criterion. 
The neutrino mass matrix in the {i'i,N'^, S'^) basis can be written as 



^ mjj m^^ 
ruD Mn 
\ms Mr fi J 



(2.1) 



In the literature many variants of this model have been considered. 
Inverse Seesaw 
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The conventional inverse seesaw models assume the terms mg and in Eq. 2A to be 
zero. Then the mass matrix takes the form, 

/ m'^ \ 
M^= mo Ml . (2.2) 
V Mr fl J 

The model is lepton number conserving in the limit ^ tending to zero. The minimal inverse 
seesaw model considered in the literature [39] consists of Svl + "^^R + 25. The model with 
Svl + + 15 is a 5 X 5 matrix with rank 3. Thus there are two zero eigenvalues of this 
matrix which is not consistent with low energy phenomenology. The model consisting of 
Svl + 2Nji + 15 is a 6 X 6 matrix with rank 5. Thus there is one zero eigenvalue. However, 
this belongs to the {Nr, 5) block and hence this scenario is not considered if one assumes 
that there are no light singlets. In Eq. ( |2.2| ) the Majorana mass term of Nr is not consid- 
ered. In principle this can be included [40,41], although this does not change the structure 
of the effective light neutrino mass matrix at the leading order [42,43]. 
Linear Seesaw 

In the so called linear seesaw models [25-27] one retains the u — S term in the Lagrangian 
through the Yukawa coupling matrix Yg and makes the ^ and the M^v term to be zero. In 
these models lepton number violation stems from the term containing Yg. The correspond- 
ing mass matrix is, 



/ 



M^ 



T 

mi) rrig 



niD Mr\ . (2.3) 
\ms Ml 

In the limit Mr » mD,iTis the above mass matrix can be diagonalized using the seesaw 
approximation and in the leading order the effective light neutrino mass matrix mught can 
be expressed as 

might = m^MR'^ms + rvLsMR'^mD- (2.4) 

Since this contains only one power of the Dirac mass term it is called linear seesaw. 

One can make an order of magnitude estimate of the various terms to check the con- 
ditions required to get rriu ~ 0.1 eV. Assuming typical values niD ^ 100 GeV (Yukawa 
coupling strength Y^ ~ 0(1), v ~ 100 GeV) and Mr = 1 TeV one needs Ys ~ 10"^^ 
In the heavy sector we get two degenerate neutrinos of mass ~ TeV. The minimal model 
consists of adding just two singlet states A''^ and 5. The rank of the 3+1+1 mass matrix 
is 4 corresponding to one zero mass eigenvalue. The Majorana mass term M^r can also be 
included which would lift the degeneracy between the heavy states. However, the contri- 
bution of this term to the light neutrino mass matrix is sub-dominant [42,43]. 
Inverse + Linear Seesaw 

It is also possible to keep both the terms and /i in the Lagrangian. Then in the limit 
Mr >> m£)^ms and in the leading order the effective light neutrino mass matrix rringht 
can be expressed as 

might = -rn]^Y72^D + rn]^^ms + rns^mr). (2.5) 

M^ Mr Mr 



- 5 - 



In this case, for Mr ~ TeV, one needs ~ 10"^ GeV and Ys ~ 10"". This hybrid scenario 
aUows one to reconstruct Yy and the combination Yg — jM^^'^- Thus reconstruction of Ys 
requires another unknown parameter, ^ [32]. In our subsequent discussion we assume ji to 
be zero and consider the hnear seesaw option. 



3. Minimal Linear Seesaw Model 

The Minimal Linear Seesaw Model (MLSM) is defined by the mass matrix in Eq. (| 
with just 2 singlet states. Then the entries Mji are numbers instead of matrices and the 
dimension of the full matrix is 5 x 5. 
This can be recasted as, 

M.=( (3.1) 

where m'^ = (m^,m^) Then defining M as 

The neutrino mass matrix can be diagonalized by a 5 x 5 unitary matrix Uq as 

U^WUo = M^"'?, (3.3) 

where M'^*'^^ = diag(mi, m2, ms, Mi, M2) with mass eigenvalues nii [i = 1,2,3)) and Mj 
(j = 1,2) for light and heavy neutrinos respectively. Following standard procedure of 
two-step digonalization Uq can be expressed as [44] 

where W is the matrix which brings the full 5x5 neutrino matrix, in the block diagonal 
form 




miight 

'^heavy 



(3.5) 



Uu = diag{U, Ur) diagonalizes the mass matrices in the light and heavy sector appearing 
in the upper and lower block of the block diagonal matrix respectively. Ul in Eq.(^j^) 
corresponds to Upmns which acquires a non- unitary correction (1 — e/2). The eigenvalues 
(Ml, M2) are obtained as {—Mr, Mr) corresponding to degenerate neutrinos with opposite 
CP parities. The negative sign in the mass eigenvalues can be absorbed in the phases of 
the diagonalizing matrix Ur giving, 

1 / 1 \ 

(3.6) 
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e and e' which characterize the non-unitarity, are given by 



M-^m'^ml{Ary . (3.7) 



Since Eq.(3.1) is in the standard seesaw form it is straightforward to obtain the light 



neutrino mass matrix 

might = m'^M~^m'jj (3.8) 

in the hmit >> m£),ms. Now inserting the expression for m^, the hght neutrino 
mass matrix is the same as that in Eq. ( |2.4| ). Note that the complete mass matrix ( |2.3| ) for 
the minimal model has 7 phases out of which 5 can be rotated away by redefinition of the 
fields. Thus there are 2 independent phases in this matrix. We choose the basis in which 
Alji is real and attach the phases to the elements of Yj^ and Yg. Since Mi, for this case is 
of rank 4, there is one zero eigenvalue. Thus one of the light neutrino states is massless 
and the two remaining masses are completely determined in terms of the two mass squared 
differences measured in oscillation experiments. 

It is very interesting to note that mught for this case is determined in terms of two 
independent vectors 

Y, = yA; Ys = ysh (3.9) 

where a and b are complex vectors with unit norm, y^, and yg are the norms of the Yukawa 
matrices and Ys, respectively. This feature allows one to completely reconstruct the 
Yukawa matrices Yi, and Ys in terms of the oscillation parameters as [32] , 

• Normal Hierarchy (NH): (mi < m2 < m^) 

Ys = ^ (yrr^ ui - c^l) (3.10) 

with 



yi + r- 

UiS are the columns of the unitary matrix U that diagonalizes the light neutrino mass 
matrix (mught) above and r is the ratio of the solar and atmospheric mass squared 
differences 



© 



(3.12) 



Inverted Hierarchy (IH): (ms << m2 ~ mi) 



Ys = ^ (71+7 ul - Uf) (3.13) 
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with 

We use the following form for U 
( 



(3.14) 



Cl2 Cl3 
C23 •S12 - S23 Sl3 Cl2 6* 
\ S23 Si2 - C23 Si3 Ci2 6**^ 



C23 Cl2 - 
-S23 Cl2 



S12C13 
- S23Sl3'Sl2e 
C23 Sl3 Sl2 e 




(3.15) 



where Cij = cos 6ij, Sij = sin 9ij, 5 is the Dirac CP phase. For the Majorana phase matrix 
P, we use P = diag(e-*", e*", 1). Note that in this case since one of the mass eigenvalues 
is zero there is only one Majorana phase. The 3a ranges of the oscillation parameters are 
tabulated in Table |l] [45]. 

From the above forms of the Yukawa matri- 
ces it is evident that these are completely deter- 
mined in terms of the masses and mixing angles, 
two unknown phases and the norms of the Yukawa 
couplings yu and Us- 

The minimal Type-I seesaw model also con- 
sists of three left-handed and two gauge singlet 
right-handed neutrinos [28-31]. However both the 
right-handed neutrinos are assumed to have the 
same lepton number. Thus both have lepton num- 
ber conserving Dirac type coupling with the light 
state. In order to have the right-handed neutrinos 
of this model at TeV scale one needs to have small 
values of the Dirac coupling Y^j ~ 10"^ unless one 



Parameters 


3a range 


Am| [10-5 eV^] 


7.12 - 8.20 


Aml,^ [10-3 eV2] 


2.31 - 2.74 
2.21 - 2.64 


sin^ 6*12 


0.27 - 0.37 


sin^ 6*23 


0.36 - 0.68 
0.37 - 0.67 


sin^ 6^13 


0.017 - 0.033 


6 


- 27r 



Table 1: Present 3a range of neutrino 
oscillation parameters. The upper (lower) 
row corresponds to normal (inverted) hi- 
erarchy. Values of Amj^ and sin^ 612 are 
hierarchy independent [45] . 



allows for fine tuning leading to m/jMr 



m 



D 







[11-13]. Since this coupling is not lepton number 
violating, its smallness cannot be explained natu- 
rally. Also for such small values of the coupling, collider signals would be suppressed even 
though the mass of the right-handed neutrino is at TeV scale. 

4. Vacuum stability of the Higgs potential 
4.1 Higgs mass and Vacuum stability in SM 

The tree-level potential of the Higgs field in the the Standard Model(SM) is given as 

y($) = A (^^V) ^ - mVV- (4.1) 

This receives quantum corrections from higher order loop diagrams. As a result the Higgs 
quartic coupling (A) runs with energy scale. The Renormalization Group (RG) equation 
for the Higgs quartic coupling A can be expressed in general as 



dX_ 
^df, ^(167r2)- 



(4.2) 
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where i denotes the z*'' loop. Assuming SM to be vahd upto Planck scale the /3 function 
calculated upto 1 loop is given as, 



/3« = 24A^ - i^-gf + 9gi)X + ^gf + ^gfgl + + 4TA - 2Y, (4.3) 



where, 



T = Tr 
y = Tr 



SYJYu + SYd^Yd + V,] , (4.4) 
S{YjY^f + 3(yrfty^)2 + (y^ty^)2l (4,5) 



In the above equations, gi denote the gauge coupling constants with i = 1, 2, 3 correspond- 
ing to U{1), SU{2) and SU{3) groups respectively. The above equations include the Grand 
Unified Theory (GUT) modified coupling for the U(l) group. Yy with f = u,d,l represent 
the Yukawa coupling matrices for the up and down quarks and the charged leptons. The 
running behavior is controlled mainly by the top quark mass rrit which drives A towards 
more negative values in the low Higgs mass region. The running of the top Yukawa is 
governed by the following equations 



^(1) _ Y 



(4.6) 



In the numerical work we have used two loop Renormalization Group Equations (RGE) 
for A, the Yukawa and the gauge couplings [46-50]. As discussed earlier the constraints 
from vacuum stability and perturbativity limits Higgs mass in the range 126-171 GeV. 
Therefore if the scalar particle observed by the CMS [1] and ATLAS [2] collaboration is 
assumed to be the Higgs Boson then the reported mass is near the lower bound obtained 
from vacuum stability condition. In the upper panels of Fig. |l], we plot the running of A 
as a function of the renormalization scale for illustrative values of Higgs mass {irih), top 
mass {rrit) and strong coupling constant (a^ = g'^/Air). The allowed range of values of rrit 
(173.2 ± 0.9 GeV) is taken from [51] and that of (0.1184 ± 0.0007) is taken from [52]. 
The Higgs mass has been varied between 124.6 — 126.6 GeV, combining the ATLAS and 
CMS results, Eq. 

We have included the corrections to incorporate the mismatch between the top pole 
mass and MS renormalized coupling. This is given as [36], 

yt{mt) = ^^{l + St{mt)) (4.7) 

V 

5t{mt) denotes the matching correction at top pole mass. We include the QCD correc- 
tions upto three loops [53] while electroweak corrections are taken upto 2 loops [54]. We 
also include suitable matching conditions for MS renormalized A and the Higgs mass at 
fi = mt [55]. The threshold effect due to the top mass is included. The plots demonstrate 

(i) Higher values of top mass drives A towards more negative values at a lower renormal- 
ization scale. 

(ii) The uncertainty in the value of the strong coupling constant also affects the running 
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ot "0 1191 
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a, = 0.1184 




m, = 173.2 (GeV) 




a, = 0.1177 




m, = 174.1 (GeV) 





10^ 10^ 



H(GeV) 



mi, = 126.6 (GeV) 
yv = 0.15 
a, = 0.1191 
m, = 172.3 (GeV) 

a^ = 0.1184 

m, = 173.2 (GeV) 

a^ = 0.1177 

m, = 174.1 (GeV) 



H(GeV) 



Figure 1: The upper panels show the variation of A with the renormahzation scale for fixed values 
of the parameters (m/i, as, rrit). The lower panels show the changes after including the Dirac Yukawa 
coupling parameter y^. 



of A in a direction opposite to the top mass effect. Higher values of Os aids in keeping A 
positive. 

The plots corroborate the fact that for lower values of Higgs mass in the range reported by 
ATLAS and CMS, the stability of the vacuum till the Planck scale may be retained only 
for certain choices of the parameters (771^,0^). For instance for nih = 124.6 GeV, A = 
is achieved around 10^*^ — 10^^ GeV for higher values of mt and lower values of a^. A 
non- negative A till the Planck Scale is achieved only for the mass of the top (a^) near its 
lowest (highest) range. For a Higgs mass of 126.6 GeV vacuum stability condition at Mpi 
can be attained for a wider range of values of the other parameters. 



4.2 Vacuum Stability in the Minimal Linear Seesaw Model 

The MLSM contains two degenerate singlet neutrinos at TeV scale. Presence of these fields 
modify the SM Renormahzation Group Equations (RGEs) for the Yukawa couplings and 
the Higgs self coupling, for energies higher than the mass of the singlets. Including the 
corrections due to the neutrino Yukawa couplings upto one loop, the modified (3 function 
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governing the running of A is given as, 



+ 4TV(y/ty^')A _ 2Tv[{Y,'^Yl) 



(4. 



Where Yl'^ = [YJ, Ys)- "^^^ ^'^^ loop /3 functions corresponding to the Yukawa couphngs 
Yu, Yd and 1/ also acquire additional factors containing Yj'^Yl [56]. Finally, one needs to 
include the RG running of the coupling Yl which is governed by the following equation: 



2 (i^l 



Y' 



-Y 'V 



Yi^Yi + T- 



20 



9l 



9 2 
792 



In this case the quantity T is defined as, 



Tr 



3YJYu + SYd^Yd + Yi^Yi + y/ty^ 



(4.9) 



(4.10) 



RG equation for neutrino Yukawa coupling is taken upto one loop [56]. 

The Yl dependence of the beta function of A is in terms of Tr [Y^'^Y^] and TV [{Y^'^Yl)^] 
only. From the parameterization of Yi, and Ys, we find 



Tr 



Y 'V 



2 , 2 
Vu + Vs- 



Tr 



Y 'Wy 'W 



yt + '^yWsP^ + yt^yt, 



(4.11) 



since Us «yu- 

Remarkably these quantities do not depend on the neutrino oscillation parameters. In 
addition under the approximation of ys « y,^, there is no dependence on mass hierarchy 
as well. In Fig. |l| we show the effect of inclusion of this term on the running of A. As 
expected, A becomes negative at a lower renormalization scale in presence of the seesaw 
term. 

It is to be noted that the vacuum stability condition, X{Mpi) = 0, used above cor- 
responds to the tree level potential. Considering the one loop effective potential for the 
Higgs field, the stability of the electroweak vacuum demands A = at Mpi, where A is the 
one loop corrected self coupling* defined as [57, 58] 



A = A 



1 



3, 2 2^2 A , [al+al] 



Vt 



1 



(4.12) 



In Fig. I we give the plot of the allowed region of function of the Higgs mass 

for fixed values of top mass and the strong coupling constant. The two curves with the 
same line type(color) that are close to each other correspond to the condition \{Mpi) = 
(the upper curve) and \{Mpi) = (the lower curve). The region below the curves satisfy 
the condition of vacuum stability, yu = corresponds to the SM. The first panel is for 
rrit = 173.2 GeV. For this value of top mass, only a small allowed region in y^ is obtained 



'Contribution of heavy neutrinos towards the efTective potential has not been considered. 




1 25 1 25.5 126 1 26.5 1 25 1 25.5 126 1 26.5 



mf,(GeV) mt,(GeV) 

Figure 2: The allowed region of yi, with varying Higgs mass. The region below the curves is 
allowed. Among the two same type (color) lines, the upper line corresponds to A(Mp;) > and the 
lower line corresponds to X{Mpi) > 0. 




172.4 172.6 172.8 1 73 173.2 173.4 172.4 172.6 1 72.8 173 173.2 173.4 

m((GeV) m((GeV) 



Figure 3: Same as in Fig. ^ with varying top mass. 

for Higgs mass > 126.5 GeV and ag = 0.1191, considering the condition X(Mpi) = 0. From 
this we obtain an upper bound yu ^ 0.1. If we ignore the one-loop correction to the self 
coupling then the allowed region is slightly larger. The 2'^'^ panel corresponds to a lower 
value of rrit = 172.3 GeV. In this case we obtain the upper bound Ui, ^ 0.23 — 0.3 depending 
on the value of Og. This bound is obtained for a Higgs mass of 126.6 GeV. For lower values 
of Higgs mass, the allowed yu values are correspondingly lower. 

Fig. ^ shows the variation of yi, with mt for fixed values of Higgs mass. Fig. ^ show 
the variation of yi, with a^. In general, larger allowed regions are obtained for higher values 
of Higgs mass, lower values of top mass and higher values of a^. From the above figures 
we can have an overall upper bound on the value of yi, as 

y, < 0.3. (4.13) 

The above plots are obtained by keeping M/j fixed at 1 TeV. In Fig. |5| we show how the 
upper bound on yi, depends on the scale of Mr. 
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Figure 5: The left-panel shows the upper bound on yi, as a function of the right-handed neutrino 
mass from consideration of vacuum stability. The upper set of lines are for mt ~ 172.3 GeV and 
the lower lines are for nit = 173.2 GeV. The red (green) solid (dashed) line in the right panel shows 
the upper limit on j/g (same as the lower limit of y^) obtained from the measured value of \/Am^^ 
for NH (IH). Also shown are the lower bound on ys from vacuum stability and neutrino masses. 
See text for details. 

We see that variation of Mji within a few TeV (which is our range of current interest) 
would not change the bound on r/i, drastically. In fact this trend continues even if Mr 
is increased to higher values. At Mr = 10^^ GeV, the upper bound on y,^ obtained is, 
Vu ~ 0.41. Note that for higher values of Mr, ys needs to be increased to keep the neutrino 
mass in the desired range. Beyond 10^^ GeV the contribution from the ys term starts 
getting significant and hence this needs to be included in the RG evolution. 

5. Constraints from Neutrino Mass 

In the linear seesaw model lepton number breaking is introduced by the term ms while 
the term mo conserves lepton number. Since the absence of the term ms enhances the 
symmetry of the Lagrangian it is natural to assume that the lepton number violating term 
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is much smaller than the lepton number conserving term. This consideration puts an 
absolute upper bound on ys which is 



{ys)max < Vu- (5.1) 

This also defines a lower limit on y^^ since for smaller values of t/s, yi, has to be higher so 
that the relation y^ysv'^ /Mji ~ is satisfied. Thus the above equation can be made to 
imply a stronger statement which is 

{Us^max ^ {yu)min' (5-2) 

On the other hand, the upper bound on y^ at each M/j, as obtained from the vacuum 
stability condition can be used to obtain a lower bound on , again from the standpoint of 
producing correct neutrino masses. From these considerations it is possible to set a lower 
bound on y^, and lower and upper bounds on for both NH and IH. The right panel in Fig. 
m displays these bounds on ys as a function of M/j incorporating the constraints on neutrino 
mass squared differences from oscillation data for NH and IH. The mass eigenvalues can 
be expressed as, 

mi = 0, ^2 = "^{1-9). -3 = ^(1 + P) {NH), (5.3) 
where m2 and are related to the mass squared differences (Am|- = m?- — mf) as. 



1712 = ^J Amj-^ = ^ Ami, "^3 = ^ Am|2 + Am^^ = ^Am^^^ + Am|. (5.4) 
This gives us the following equations for NH: 



(1 - p) (2^) 2Mr ■ 
Taking p ~ 1 — for NH, the above equations give a single equation: 



(5.5) 



y.Vs « MR^Aml^Jv\ (5.6) 

Using Eq. ( [5^ ) in Eq. ( |5^ one gets (with {ys)max = {yu)min), 

{ys)max « ^/M^{Amlt^f'^/v. (5.7) 

This defines the upper bound on (or the lower bound on y^) for each Mr. The solid 
(red) line in the right panel in Fig. ^ displays these bounds as a function of Mr. The 



plot is obtained using the exact Eq. (5^) and the uncertainty in the values of Ami, and 



Am^j^ gives a thin band which is displayed in the inset. On the other hand since we know 
{yu)max at each Mr from consideration of vacuum stability, the above equation gives us a 
lower bound on y^ as 



Mr^/A 



y.>^|f— (5.8) 

^ yyujmax 



2 
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This is also shown in Fig. |5|. 

For IH the mass eigenvalues are given as, 

The absolute masses can be expressed in terms of mass squared differences as, 



mi = y Amfg = ^J Am^^^, ma = ^ Am| + Am^^^, (5.10) 

and we get, 



(1 - p) (2^) 2Mr 
Taking p r/4: for IH, we get a single equation 



(5.11) 



y,ys^2MR^/\ml,Jv\ (5.12) 

Again this gives an upper bound on which is the same as the lower bound on y^. This 
is shown in the right panel of Fig. ^ by the dashed (green) line. Similarly one can also use 
this equation and the upper bound on y^ obtained from vacuum stability condition to put 
a lower bound on yg which is displayed in Fig. ^ 
The bounds on for Mr = 1 TeV : 

(3.3 - 3.6) X 10"^^ <ys< (1-0 - 1.04) x 10"^ NH (5.13) 
(5.5-6.0) X 10"^^ <y, < (1.28 - 1.34) X 10"^ IH (5.14) 

Thus using constraints from vacuum stability and neutrino masses, one can determine the 
range of the unknown coupling parameters y^ and y^. 

6. Non-unitarity of Neutrino Mixing Matrix and Lepton Flavor Violation 

The mechanism of neutrino oscillation has already indicated that there is flavor violation in 
the lepton sector. The question arises if there can be flavor violating decays in the charged 
lepton sector. In typical Seesaw models this rate is very small because of the smallness of 
the light-heavy mixing. However in TeV scale seesaw models, since lepton number violation 
(LNV) is separated from the scale of lepton flavor violation (LFV), this may not be the 
case [59]. It is well known that LNV is due to the dimension 5 operator whereas the LFV 
can be related to the dimension 6 operator. In general, the flavor structure of the coupling 
strengths of these operators are not correlated. However, in MLSM such a relation can be 
established from the hypothesis of minimal flavor violation [60]. In this model, the mass 
and gauge eigenstates are related as, 

{vL v'') = U {vl' v^') (6.1) 

where, 

= (n% 5^) (6.2) 
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The leptonic part of the charged current interaction in the gauge basis is given as , 

9 

This can be expressed in the mass basis as, 

3 2 



^CC = ^ ^aL7^iyaL W + h.C. 



(6.3) 



-cc 



_g_ 

V2 



a=e,fi,T i=l j=l 

The PMNS matrix is defined as 



E T.T.[^'»L^^^^li(^L)a^ly'u + {VUNf,^}W>']+h.c. (6.4) 



f/pMNS = U1 (l-^e]U 



(6.5) 



where Ui is the unitary matrix which takes the left-handed charged lepton fields to their 
mass basis and other quantities are defined earlier. As we are working in a basis, where 
charged lepton mass matrix is diagonal, Ui is being taken as unity. We note that Upmns 
is non-unitary and the correction to unitarity is proportional to e/2. 

In this section we consider the branch- 
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> 67) 
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Br(r - 


> 67) 
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Br(r - 




< 4.4 X 10-*^ 


Br(;U - 


> 3e) 


< 1.0 X 10~^2 


Br(r - 


> 36) 


< 2.7 X 10^*^ 


Br(r - 




< 2.1 X 10"^ 


Br(T ^ 




< 1.7 X 10-^ 


Br(r ^ 


ee/i) 


< 1.5 X 10"^ 



ing ratios of LFV decays in MLSM. In 
view of the recent measurement of ^13 
the branching ratios now can be stud- 
ied in terms of the CP phases. In addi- 
tion, from the experimental upper lim- 
its on LFV processes one can obtain 
constraints on the parameter y^/Mji as 
a function of the CP phases. When 
combined with the upper bound on y^, 
from vacuum stability as a function of 
Mji, the parameter space can be fur- 
ther constrained. Table |2| lists the ex- 
perimental constraints coming from the charged lepton flavor violating decays [61]. 

Nevertheless, in this section we will concentrate only on the constraints coming from 
Br(^ — )• 67) since this is the most constraining as can be seen from Table |2|^. 
Branching ratio for the process, ;U — )• 67 is given by [62] 



Table 2: Various experimental constraints from 
charged lepton flavor violating decays [61]. 



Br (// 



3a 



yeiVlfix) 



(6.6) 



where, x 



w 



m 



and 



fix) 



X (1 — 6x + 3x'^ + 2x'^ — Ox^lnx) 
2(1 -x)4 



(6.7) 



^Although the — ^ 3e conversion in nuclei has a better constraint, the sensitivities are lowered by the 
uncertainties in nuclear form factors [21]. 
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f{x) is a slowly varying function of x ranging from to 1 for x between to infinity. The 
light-heavy mixing matrix V is define in Eq. 3.4. The current experimental constraint on 
this is [63] (see, table ^ 

Br (^ ^ 67) < 2.4 X 10-^2^ (6.8) 

Using the parameterization of and Ys in Eq. (|3.1C ) and (3.13), we obtain, for normal 
hierarchy 

3a y'^v'^ 



Br (fi — 67) 



Svr 4M^ 



r(^) 



2,.(i/4) 



+ 0{ys,{V^,si3f) . (6.9) 
In the above expressions and in subsequent part, we have used the following notations 



sin0. 



si: 



= cos9- 
c4i 



= cos4:6ij, 



= cos (a + 6), s2ij 
c2a = cos2q, c4o 



sin20,- 



c2i 



cos4q; etc. 



cos26'. 



.10) 



The above equation in conjunction to the upper bound on the i?r(/i — )• 67) can be 
used to put an upper bound on y^/M^ as, 

1/4 



Vu/Mr < 



2.4 X 10 



-10 



3av^Pix)G^Hir,9ij,a,6) 



(6.11) 



Here, the factor {r, 9ij,a, S) contains the oscillation parameters. The upper bound 
on yiy/M^ varies in a range depending on the values of the CP phases. The minimum value 
of the upper bound occurs at a + 5 = while the maximum occurs at a + 6 = n. This 
is reflected in the top panels in Fig. |^ where we display the allowed values of y^/Mn as a 
function of the CP phases'''. The left most panel displays the variation of the upper bound in 
yyjMji as a function of the Majorana phase a. The solid(red) line corresponds to the Dirac 
phase = while the dashed (green) line is for 5 = vr. The other oscillation parameters 
are marginalized over the 3(T range in Table 1 to give the maximum and minimum value 
of the upper bound on yy/M^. For other values of oscillation parameters the upper bound 
would lie somewhere between these two values. From the figure it can be inferred that the 
maximum value of the upper bound on y^/Mji is 

< 0.00049(GeV-^), (6.12) 

which occurs for 5 = and a = tt for NH. 

For IH, the branching ratio can be expressed as, 

-2si3 (2s5s2„(l + c2„s2i2) + c2„C5(2c2i2 + c2„s4i2)) S223] 
+0{ys,{V^,si3f) (6.13) 



■''In this plot we have taken f{x) to be unity. For M_r = 200 GeV, there will be a multiplicative factor 
~ 1.3. As Mr increases, this factor tends to become unity. In Fig. ^, we have included the exact value of 
f{x) at each Mr. 
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Figure 6: The upper panels show the allowed regions of Ui^/Mn vs the CP phases a and 5 for NH 
while the lower panels are for IH. The area below each curve is consistent with the experimental 
upper bound on the rate oi — >■ 67. The two lines of the same type (color) correspond to the 
maximum and minimum value of the upper bound obtained by varying the oscillation parameters 
over their 3cr range. 



From this one can again get an upper bound on y^/Mfi. This is displayed in the lower 
panels in Fig. ^. As for NH the value of the upper bound depends on the CP phases. The 
maximum allowed value in this case is 

2/,/Mfi<0.7(GeV-i), (6.14) 

which occurs for 5 = 0, q = as can be seen from the figure. The lines of same line type 
(color) corresponds to the upper bounds including the uncertainties in the masses and 
mixing parameters. 

The maximum value of yu/Mn as obtained above can be used to retrieve the maximum 
value of yi, for each Mr. This is shown in Fig. |7[ Note that while extracting the bound 
on yiy for a particular Mr from the figure one has to be careful to ensure that the per- 
turbativity bound on y^ 1) is not violated. We also superimpose the bounds obtained 
from consideration of vacuum stability in this figure. The area to the left of the shaded 
bands is disallowed from the constraint on the branching ratio fi — t- ej. These bands are 
obtained for fixed values of the CP phases {6, a). The band for each combination of CP 
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Figure 7: The allowed regions of j/™"^ as a function of Mji from the combined constraints of 
Br(i.i — ^ 67) and vacuum stability. The area to the right of the curved lines are allowed from 
experimental bound on i?r(/i — >■ 67) while the area below the slanting lines are allowed from the 
constraint on vacuum stability. The vacuum stability bound shown here correspond to mt = 172.3 
GeV and = 0.1191. The two (very) close line corresponding to nih — 126.6 GeV is due to 
X{Mpi) = (upper one) and X{AIpi) = (lower one). For rrih — 124.6 GeV the hue corresponds to 
X{Mpi) = 0. No allowed region is obtained in this case for X{Mpi) = 0. 



phase is obtained by varying the oscillation parameters in their current 3a range. The area 
below the slanting lines are allowed from the constraint on vacuum stability. The figure 
shows that the constraints from Br(/u — )> ej) can sometimes further constrain the value of 
yiy as obtained from vacuum stability. For instance for Mr = 200 GeV and NH the con- 
straint from Br(/i — t- cy) restricts yu to be < 0.13 for values of CP phases {S,a) = (0, vr). 
For {6, a) = (0, 0) the maximum allowed value of is lower. For other combinations of 
CP phases the bands lie anywhere inside or between the two shaded regions. Thus if we 
consider all possible values of CP phases then only the region marked disallowed is not 
compatible with the constraints from Br(^ — ?■ 67) for NH though it was consistent with 
vacuum stability constraints. 

For IH and {5, a) = (0, 0), the hatched region extends all the way upto Mr = 200 GeV 
and there is no significant constraint from /i — )■ 67 given the present uncertainty on the 
neutrino oscillation parameters. However for the green hatched region corresponding to 
{5, a) = (0,-7r/4) the region to its left is disfavored and y^"^ is constrained to lower values 
as compared to the bound from vacuum stability. But if we consider all possible values 
of CP phases then this area becomes allowed. Hence we conclude that, given the present 
uncertainty of oscillation parameters and the CP phases, no improvement over the vacuum 
stability bound can be obtained on {yu)max from Br(/i — )• 67) for IH. 

We note in passing that in this type of models the Higgs boson can decay to two 
neutrinos of which one is heavy and the other one is a light neutrino, as long as the heavy 
neutrino is lighter than the Higgs boson. This has been studied in the context of inverse 
seesaw models and put constraints on the Yukawa coupling y^, to be y^ ^ 0.02 for Mr ^ 
120 GeV from the experimental data on the channel h — )• WW* — t- iivi' [16,68]. For larger 
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masses of the heavy neutrino current Higgs searches do not provide any constraint on the 
parameter space. 

On the other hand, the search for heavy singlet neutrinos at LEP by the L3 collab- 
oration in the decay channel N — t- eW showed no evidence of such a singlet neutrino in 
the mass range between 80 GeV {\Vai\'^ < 2 x IQ-^) and 205 GeV {\Vai\'^ < 1) [64]. V^i 
is the mixing parameter between the heavy and light neutrino. Heavy singlet neutrinos 
in the mass range from 3 GeV up to the Z-boson mass (mz) has also been excluded by 
LEP experiments from Z-boson decay upto |V^i|^ ~ 10~^ [65-67]. In the light of these 
experimental observations we have chosen the parameter Mr to be greater than or equal 
to 200 GeV in this study. 



7. OZ//3/3 decay 



The half life for neutrino-less double beta decay in presence of heavy singlets is given 
by [69,70]. 

2 



-(1/2) 



G 



mi 



Ul . m,+ < p2 > Ik 



where < > is given by 



< p > -- 



-rriemp 



M 



N 



(7.1) 



(7.2) 



Ail, and Ai n denote the nuclear matrix elements corresponding to light and heavy neutrino 
exchange respectively. The values of the parameters are taken as [69] G = 7.93 x 10~^^ 



yr-i, < p2 >= _(i82 MeV)^. 



The first term in Eq. (|7. 1| ) is the usual contribution from the left-handed neutrinos. 
The second term denotes the contribution of the singlet neutrinos. The matrix V is defined 
in Eq. ( |3.4] ). Taking the most general form of the matrix as 



nidi rudz 
, nisi ms3 , 



and Ur and M as defined in Eq. (|3.6D and (|3.2|) respectively, we obtain 



el 



-(m*i -m2i),V; 



1 



Then the contribution from the heavy part is 2 < > m*^m^j^/M^ ~ 10~^mj. Thus this 
contribution is negligible as compared to the contribution from the light sector which is 
~ rrii- 

Therefore, Oz^/3/3 is due to the light neutrinos only and the effective mass is defined as 



-{1^*81 +m*di)- 



(7.3) 



(7.4) 



(7.5) 



Since in this case the lightest mass is zero one can plot the conventional plots of effective 
mass as a function of the unknown CP phases for both hierarchies. 
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□ Inverted Hierarchy 
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I ~1 Inverted Hierarchy 



Figure 8: The effective mass governing 01/(3/3 as a function of tlie Majorana pliasc a for NH (dark 
(red) shaded curve) and IH (light (green) shaded curve). The left panel is ioi 6 — while the right 
panel is for (5 = 7r/2. 



For NH the effective mass ruee in the limit of the smallest mass mi — )• is given as, 

(7.6) 



\mee\NH = \ Am, 



^2 

"atm 



V-rslA.e'-'^ + she-''' 



The maximum is obtained for {a, 6) = (0,0) or (7r/2,7r/2) while the minimum occurs for 
{a, 6) = (0, 7r/2) or (7r/2,0). This is reflected in Fig. ^ by the dark (red) shaded curve 
which represents the effective mass governing 01/(3(3 as a function of the Majorana phase a. 
The shaded portion is due to the 3a uncertainty in the oscillation parameters that appear 
in the expression of effective mass. The left panel is for 5 = and the right panel is for 
6 = 7r/2. The cancellation condition is 

^/r sin'912 = tan^ 6*13, (7.7) 

which is not satisfied for the current 2>a ranges of parameters and therefore the effective 
mass does not vanish which is also seen from the figure. For IH the smallest mass is ms 
which is zero in this model and the effective mass is 



<tm{42che-'''' + siA,en- (7.8) 

For IH the effective mass is independent of the Dirac phase 5. The maximum of |mee| 
occurs for a = 0, 7r/2,7r and the corresponding expression is. 



'13 y ''atm' 

The minimum value is obtained for a = 7r/4, 37r/4 as. 



^mltm- (7.9) 



\<"\min = 4cos20i2^Am2,^. (7.10) 

This is seen from Fig. ^ by the light (green) shaded curve, rn-ee for IH is in the range 
accessible to future neutrinoless double beta decay experiments. 
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8. Collider Signatures 



As mentioned earlier, if the heavy singlet neutrinos have mass less than the Higgs boson 
mass, then the Higgs boson can have new decay modes [16]. For example, the Higgs boson 
can decay into h vN . Now, the singlet neutrinos can decay into IW and uZ through 
the mixing between the heavy neutrinos and the light active neutrinos. At the LHC this 
will lead to final states such as pp ^ h ^ i~^£~ + EJt , where £ = e, /i. Note that these 
final states will depend on the Yukawa couplings and one can put bounds on these Yukawa 
couplings from the existing LHC data on these types of final states [16]. 

We have considered the singlet neutrino to be heavier than the Higgs boson. In this 
case one has to look at the 3-body decay modes of the Higgs boson through the virtual 
heavy neutrino to have similar final states. Obviously, in this case the constraints on the 
Yukawa couplings will be much less restrictive. In our model we have obtained upper bound 
on the Yukawa couplings y^, from the vacuum stability condition and this can be used to 
test our model at the LHC by looking at the dilepton plus missing Et final states. 

One can also have trilepton plus missing Et final states at the LHC from the production 
of these heavy neutrinos [17]. For example, at the LHC the heavy neutrinos can be produced 
through the s-channcl exchange: ud — > i'^N or ud — t- i^S. N oi S can again decay 
into IW and vZ through f — N or v — S mixings. This will lead to trilepton plus missing 
Et final states at the LHC. Now, the trilepton plus EJt signal is a very clean signal for 
looking at physics beyond the standard model. In our model, the trilepton final states 
depend once again on the Yukawa couplings of our model. Using the upper bound on y^j 
obtained from vacuum stability condition, it would be possible to study the present model 
at the LHC through the trilepton channel. 

9. Conclusion 

In this paper we consider the phenomenology of the minimal linear seesaw model consisting 
of three left-handed neutrinos and two singlet fields. The two singlet fields have opposite 
lepton numbers. Smallncss of neutrino mass is ensured in this model by the tiny lepton 
number violating coupling {yg) of one of the singlets with the left-handed neutrinos. Thus, 
the masses {Mr) of the heavy singlet neutrinos can be at the TeV scale even with the Dirac 
type coupling {y^) between the other singlet and the heavy state of C(l). This permits 
appreciable light-heavy mixing in the model which can have interesting phenomenological 
consequences. The model predicts one massless neutrino and hence there is only one 
Majorana phase. The great advantage of this model is that the Yukawa matrices can be 
fully reconstructed in terms of the oscillation parameters apart from the overall coupling 
strengths y^ and ys- We show that consideration of the vacuum stability of the scalar 
potential can constrain the coupling strength y^ as y^ ^ 0.3 for 200 GeV ^ Mr ^ 1 TeV. 
However, it should be noted that this bound depends on the value of the strong coupling 
constant (a^), the top quark mass (m^) and the Higgs boson mass (mii). Hence, this upper 
bound can be somewhat lowered depending on the experimental uncertainties on these 
quantities. The assumption of tiny lepton number violation implies (7/5)^^,02; ^ {yv)min- 

This 
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allows one to obtain an upper (lower) bound on (y^) in conjunction with the constraints 
on neutrino masses from oscillation data. On the other hand the vacuum stability bound 
on y^, together with the measured values of neutrino mass squared differences gives a 
lower bound on yg. Thus both the unknown overall coupling strengths of this model get 
constrained. 

We also obtain the non-unitary correction to the PMNS mixing matrix. Bounds on 
yu/Mji can be obtained as a function of the CP phases a and 5 from experimental con- 
straints on lepton flavor violating processes. Combined constraints from vacuum stability 
and the lepton flavor violating process — )■ 67 rule out a significant portion of the pa- 
rameter space in the (y^-Mn) plane for NH for masses of Mr ^ 600 GeV. On the other 
hand, contribution of the singlet neutrinos to the neutrinoless double beta decay process 
is insignificant. The model predicts interesting signatures at the LHC and can be tested 
using the present and future data. However, a complete collider study merits a separate 
analysis. 
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